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■ ABSTRACT 

We report on multiwavelength observations of the blazar PKS 0537—441 (z = 0.896) 
obtained from microwaves through 7 rays by SMA, REM, ATOM, Swift and Fermi mostly 
during 2008 August-2010 April. Strong variability has been observed in 7 rays, with two ma- 
H ■ jor flaring episodes (2009 July and 2010 March) and a harder-when-brighter behaviour, quite 

. - -' common for flat spectrum radio quasars and low-synchrotron-peaked BL Lacs, in 2010 March. 

In the same way the spectral energy distribution (SED) of the source cannot be modelled by 
a simple synchrotron self-Compton model, as opposed to many BL Lacertae objects, but the 
addition of an external Compton component of seed photons from a dust torus is needed. 
The 230 GHz light curve showed an increase simultaneous with the 7-ray one, indicating 
co-spatiality of the mm and 7-ray emission region likely at large distance from the central 
engine. The low, average, and high activity SED of the source could be fit changing only the 
electron distribution parameters, but two breaks in the electron distribution are necessary. The 
ensuing extra spectral break, located at NIR-optical frequencies, together with that in 7 rays 
seem to indicate a common origin, most likely due to an intrinsic feature in the underlying 
electron distribution. An overall correlation between the 7-ray band with the i?-band and K- 
band has been observed with no significant time lag. On the other hand, when inspecting the 
light curves on short time scales some differences are evident. In particular, flaring activity has 
been detected in NIR and optical bands with no evident 7-ray counterparts in 2009 September 
and November. Moderate variability has been observed in X-rays with no correlation between 
flux and photon index. An increase of the detected X-ray flux with no counterpart at the other 
wavelengths has been observed in 2008 October, suggesting once more a complex correlation 
between the emission at different energy bands. 

Key words: galaxies: active — galaxies: BL Lacertae objects: general — galaxies: quasars: 
general - galaxies: BL Lacertae: individual (PKS 0537—441) — gamma rays: observations 
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1 INTRODUCTION 

PKS 0537-441 (also known as HB89 0537-441, 
RX J0538.8— 4405) was selected as a quasi-stellar sour ce in 
the Parkes 2700 MHz r adio survey ( Peterson et alj Il976h and 
identified as a blazar by ISurbidge fc Hewitt C 19921) . According 
to its spectral energy distribution the source was classified as a 
BL Lac object, cate g orized as a low-energy-peaked BL Lac by 
IPadovani & Giommil ( Il995h . although the determination of the 
redshift of this source is based on broad emission lines and the ob- 
ject h a s been listed as a highl y polarized quasar dLedden fc Odelll 
Il985|; [ Sambruna et alJ Il994l) . having radio polarization up to 
3% teomesaroff et al.' 1984) and optical polarization up to 18% 
l^pev &Tapia 1988, 1990). PKS 0537-441 has also been sug- 
gested to be a transition object between classical BL La c objects 



and v i olently variable and h ighly polarized qu a sars (ICristiani 



19851; iMaraschietalJ Il985l: IXreves et all 1 19931 : loiommi et a" 



1 9951 ; iGhisellini et al J20 1 C 



Two broad emission lines at 3617A and5304A were observed 
in the optical spectrum of PKS 0537—441, and proposed to cor- 
respond to the C 111] A1909 and Mgll A2798 lines, placing the 
object at a redshif t z = 0.894 (Peterson et al. 1976; Wilkes et al. 
ll983l : IStickel et al.|[l99'3l) . Subsequently, Lewis fc Ibat j ( 120001) re- 
ported a redshift of z = 0.892 ± 0.001, based on the Balmer line 
series and the [O 111] emission line. PKS 0537—441 was also ob- 
served by the Faint Object Spectrograph (FOS) of the Hubble Space 
Telescope (HST) in 1993 July and September and, on the basis of 
their measured central wavelengths, the e mission line s detec ted in 
the HST FOS spectra were identified by IPian et al.l ( |2005|) with 
Ly Q, Si IV and C IV at z = 0.896 ± 0. 001, a value co r isisten t 
with the original redshift determination by IPeterson et al.l jl976l) . 
Based on a claimed detection of a galaxy along the line of sight to 
PKS 0537—441, it was discussed as a case of grav itational tens- 
ing jStickel. Fred. Kuhilll988l : iLewis fc Ibat3 l2000l) . in which a 
foreground galaxy lenses the radiation from the quasar, amplify- 
ing the continuum with respe ct to the line emission, but furthe r 
observations repo rted in Falomo et al. 1 ( Il992h : I Plan et al. 1 (l2002h : 
lHeidtetalJ ( l2003l) have not confirmed the presence of a foreground 
galaxy. 

Given its apparent brightness, PKS 0537—441 is an excellent 
target for studying the properties of blazars and therefore it has 
been the subject of several monitorin g campaigns from radio t o 
optical, showing remarkable variability. IPeterson & Boltoi] jl972h . 
who found Quasi Stellar Objects (QSOs) among the Parkes radio 
sources, discovered the optical counterpart of the radio source PKS 
0537—441 to be 2 mag brighter than on the Palomar O bservatory 
Sky Survey. Subsequent observations by lEggenl ( Il973[) showed a 
variation between magnitude 16.5 and 13.7 over a few months. 
The study of the Harvard Observatory photographic plate collec- 
tion iLillerlll974h over almost a century revealed long-term vari- 
ations of ~5 mag with fluctuations of ~2 mag over less than 2 
months. In addition this source displayed intraday variability at ra- 
dio ( Romero et al. 1995), and optical frequencies (Heidt fc Wagner 
ll996l : lRomero et al.ll2 002). Although PK S 0537-441 is one of the 
brightest radio selected BL La c objects (iStickel et al.lll991ll it is a 
relatively weak X-ray source dTanzi et ak 19861 : IWorral fc WilkesI 
Il990h . even if the comparison of the ROSAT X-ray spectrum with 
the Einstein and EXOSAT observations indicates sub stantial vari- 
abilit y both in intensity and spectral slope at X-rays dTreves et al.l 
1 19931 and the reference therein). BeppoSAX in November 1998 
observed the source with a spectrum and emission state consistent 
with those measured by EXOSAT and Einstein, but with a flux at 



1 keV almost a factor of two less than that detected by ROSAT 
in 1991. Moreover, the BeppoSAX spectrum suggests that a sin- 
gle emission component dominates in the energy range 0.1-30 keV 
(f ian et al. 2002). 

P revious multifrequenc y observations from near-IR to X- 
rays dMaraschi et al.l 1 19851 : iTanzi et al] Il986l) showed quasi- 
simultaneous flares at IR, optical, and X-ray frequencies. Its con- 
tinuum variations throughout the entire elec t romag netic spectrum 
have been discussed in detail bv lPianetalj d2002l) . In 2005 PKS 
0537-441 was mon itored in the optical and infrared by REM 
dOolcini et ai]|2005l) and observed by all the instruments on board 
Swift in January, July, and November (Plan et al. 2007), with a flux 
that varied by a factor of ~60 and ~4 in optical and X-rays, re- 
spectively. The l/-band and X-ray light curves measured in 2004- 
2005 were highly correlated, although with different variation am- 
plitudes. In contrast, no clear evidence of variability within a sin- 
gle night was found, even though some hints of small flares on 
day time-scales may be present. Moreover the X-ray photon in- 
dex was observed to be steady over the different epochs, as op- 
posed to the variation of the fluxes. The optical spectra collected at 
different epochs suggest the presence of thermal emission during 
the low states. Only recently, Impiombato et al. (2011) identified 
an episode of rapid variability in the J-band with a duration of ~ 
25 minutes. 

PKS 0537^41 was detected in 7 rays as 3EG J0540~4402 by 
the EG RET telescope on bo ard the Compton Gamma-Ray Obser- 
vatory dHartman etai]|l999h . The source was det ected by EGRET 
for th e first time in 1991 dMichelson et al.lll992t [Thompson et all 
Il993h and then re -observed at many successive epochs in differ- 
ent states (see also lTreves et aljl993l) . showing bright and variable 
emission in 7 rays, with a maximum flux of 9x 10^^ photons cm~^ 
s~^ during 1995 January and a peak with the temporal binnin g 
of 2 days of (20±5)x lO"'^ photons cm"^ s~^ jPian et al.ll2002h . 
More recently intense 7-ray activity from this source was observed 
by the Large Area Telescope (LAT) onboard Fermi (Tosti 200^: 
iBastieri 1 120091: Ic annonl|2010h and AGILE Gamma-Ray Imaging 
Detector dLucareUi et al.ll2010l) . PKS 0537-441 was listed in both 
the First and Second Eermi Large Are a Telescope Source catalogs 
dAbdo et al.ll2010bl : lNolan et all2012l) . 

Now the Eermi satellite, with continuous monitoring of the en- 
tire 7-ray sky and the excellent sensitivity of the LAT, gives us the 
opportunity to study the blazars in the 7-ray band not only dur- 
ing flaring states but also during low activity states, broadening 
correlated investigations of blazars over the whole electromagnetic 
spectrum. In this context the Swift satellite, with broadband cover- 
age and scheduling flexibility create a perfect synergy with Eermi 
and provide crucial data in the gap in coverage between the 7-ray 
data and the radio-to-optical data obtained from ground-based tele- 
scopes allowing a deep and complete long-term monitoring of these 
sources. 

In this paper we report the results of the LAT monitoring of 
PKS 0537^41 together with the related multiwavelength observa- 
tions across the electromagnetic spectrum. This paper is organized 
as follows. In Section 2 we present LAT data and analysis. In Sec- 
tion 3 we present the multifrequency data collected by Swift, REM, 
ATOM, and SMA. In Section 4 we discuss the 7-ray spectral and 
flux variability, while the correlation with the other energy bands is 
discussed in Section 5. The general properties of the source and a 
comparison of the multifrequency data reported in this paper with 
respect to past observations are discussed in Section 6. In Section 7 
we discuss the SED modeling for different epochs. Finally in Sec- 
tion 8 we draw our conclusions. 
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In the following we use a ACDM (concordance) cosmol- 
ogy with values given within Icr of the Wilkinson M icrowave 
Anisotropy Probe (WMAP) results llKomatsu et alj |2009^). namely 
h = 0.71, fim = 0.27, and JIa = 0.73, and a Hubble constant value 
Ho = 100 h km Mpc~^, with the corresponding luminosity dis- 
tance dL ~ 5.78 Gpc (~ 1.8 X 10^* cm). 



period. For each time bin, if the TS value for the source was less 
than 4 or the number of model predicted photons Npred < 10, a 
2 a upper limit was computed for the flux. The estimated rela- 
tive systematic unc ertainty on the 7-ray flux, which according to 
lAbdo et alj fcoiObI) reflects the relative systematic uncertainty on 
effective area, is 10% at 100 MeV, 5% at 500 MeV, and 20% at 10 
GeV. All errors reported throughout the paper are statistical only. 



2 LAT OBSERVATIONS 

The Fermi-hAT is a pair-conversion 7-ray detector, sensitive to 
photon energies from about 20 MeV to >300 GeV. It consists of 
a tracker (composed of two sections, front and back, with differ- 
ent capabilities), a calorimeter and an anticoincidence system to 
reject the charged-particle background. The LAT has a large peak 
effective area (~ 8000 cm^ for 1 GeV photons in the event class 
considered here), viewing ~ 2.4 sr of the sky with single-photon 
angular resolution (68% containment radius) of 0.6° at £ = 1 GeV 
on-axis tAtwood et al.,2009.). 

The data presented in this paper were collected in the first 20 
months of Fermi science operation, from 2008 August 4 to 2010 
April 4 (~ 600 days, from MJD 54682 to 55291). The analysis was 
performed with the standard Fermi LAT ScienceTools soft- 
ware packag43 (version v9rI5p6). Only events having the highest 
probability of being photons, belonging to the "Diffuse" class, in 
the energy range O.I~IOO GeV were used in the analysis. The in- 
strument response functions (IRFs) P6_V3_DIFFUSE were used. 
In order to avoid significant background contamination from Earth 
limb 7 rays, all events with zenith angle > 105° were excluded. In 
addition only photons detected when the spacecraft rocking angle 
was < 52° were selected. This eliminated time intervals with the 
Earth in the LAT field of view. 

The LAT light curves were computed using the maximum- 
likelihood algorithm implemented in gtlike. The 20-month light 
curve was built using 3-day time bins (Fig. [T] upper panel); for 
each time bin integrated flux, photon index, and test statistic (TS0 
values were determined. For each time bin we selected photons in- 
cluded in a region of interest (Rol) of 12° in radius, centered on the 
position of the source. In the Rol analysis the source is modelled 
with a power-law, dN/dE = (iV(r + l)E-'')/{E^+l - EI+^), 
where A'' is the normalisation, F is the photon index, E,„i,i and 
Emax are the limits of the energy interval chosen for the Likeli- 
hood analysis. 

All point sour ces listed in the First Fermi LAT catalog (IFGL; 
lAbdo et alj|2010 j) within 19° from PKS 0537-441 with TS > 50 
and fluxes above I0~* ph cm'^^s"^ were included in the Rol model 
using a power-law spectrum. The background model used to extract 
the 7-ray signal includes a Galactic diffuse emission component 
and an isotropic component. The Galactic component is parameter- 
ized by the map cube file glljem_v02.fit. The isotropic background 
component, which is the sum of residual instrumental background 
and extragalactic diffuse 7-ray background, was included in the 
ROI model using the standard model file isotropic_iem_ v02.tx{l 
In the light curve computation the photon index value was frozen 
to the value obtained from the likelihood analysis over the entire 



^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
2 The test statistic (Mattox et al. 1996) is defined as TS= -2(IogLo - 
logL), where Lq is the likelihood for the null hypothesis (i.e. no source 
exists at the given position) and Li is the alternative hypothesis (the source 
exists). 

^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 



3 THE MULTIFREQUENCY COVERAGE 

The monitoring of PKS 0537-441 by Swift ranged from optical to 
X-ray bands and was supplemented by data collected in near-IR 
and optical bands by REM and ATOM, as well as sub-mm by SMA. 
This provided excellent broadband coverage. 



3.1 Swift 

The Swift satellite iOehrels et al.l2004h performed several observa- 
tions of PKS 0537^1 between 2008 October 7 and 2010 March 
4, with all three on b oard experiments: the X-ray Telescope (XRT; 
Burrows et all |2005|) . the Ultraviolet/Optical Telescope (UVOT; 
Roming et al. 2005) a nd the coded-mask Burst Alert Telescope 
(BAT; Isarthelmv et al .' 2005). These observations were obtained 
both by means of dedicated Target of Opportunities (ToOs) and 
by activating a monthly monitoring program (PI: F. D'Ammando) 
covering the period 2009 March-2010 February. 



3.1.1 Swift/BAT 

The Burst-Alert Telescope on board the Swift satellite is a coded- 
mask telescope operating in the 15-200 keV energy range. Thanks 
to its large field of view, BAT surveys up to 80% of the sky ev- 
ery day. We selected all observations with PKS 0537—441 in the 
BAT field of view, between 2004 November and 2009 August. 
The data were processed using Heasof t package (v.6.8) and fol- 
lowing the recipes presented in lAiello et al.l J2009l) . The spectrum 
of PKS 0537—44 1 was extracted using the method presented in 
lAiello etal.l ( l2008h . 



3.1.2 Swift/XRT 

The XRT data were processed with standard procedures 
(xrtpipeline vO.I2.6), filtering, and screening criteria by us- 
ing the Heasoft package (v.6.11). The source count rate was 
low during the whole campaign (mean count rate < 0.5 counts 
s^^), so we considered only photon counting (PC) data and fur- 
ther selected XRT event grades 0-12. Pile-up correction was not 
required. Source events were extracted from a circular region with 
20-pixel radius (1 pixel ~ 2'.'36), while background events were 
extracted from an annular region centered on the source and with 
radii of 55 and 95 pixels. Ancillary response files were generated 
with xrtmkarf , and account for different extraction regions, vi- 
gnetting and PSF coiTections. 

A spectrum was extracted from each observation and fit with 
XSPEC (vl2.7.0) adopting an absorbed power law model with free 
photon index Fx and using the photoelectric absorption model 
tbabs with a neut ral hydrogen colum n fixed to its Gal actic value 
(2.9 X 10^° cm~^; lMurphv et aljl996l) . consistent with lPian ef^ 
( l2007t) . Data were rebinned to have at least 20 counts per energy 
bin to allow the minimization. 
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Figure 1. Multifrequency light curve for PKS 0537^41. The period covered is 2008 August 4 - 2010 April 4. The data sets were collected (from top to 
bottom) by Fer-mi-LAT (7 rays), Swift-XRT (0.3-10 keV), S»'i>UVOT {Wl, M2, W2, and V, B, U filters), REM {V and I bands), REM and ATOM {R 
band), REM (J, H, K bands), and SMA (230 GHz). The solid lines in the top panel represent the time intervals used in the LAT spectra extraction. 



3.1.3 Swift/UVOT 



During the Swift pointings, the UVOT jPoole et alJliooi) instru- 
ment observed PKS 0537-441 in the V, B, U, and Wl, M2 
and W2 photometric bands. The analysis was performed using the 
uvotsource tool to extract counts from a standard 5" radius 
source aperture, correct for coincidence losses, apply background 
subtraction, and calculate the source flux. The background counts 
were derived from a circular region of 10" radius in the source 
neighbourhood. The flux es were then de-redde ned using the val- 
ues of E(B-V) taken from lSchlegel et al.l (1 19981) with A(A)/E(B-V) 
ratios calculated for the UVOT filters u sing the mean interstellar 
extinction curve from lFitzpatricS (Il999l) . 



3.2 REM 

The R apid Eye Mounting (REM; IZerbi et alj|200ll ; ICovino et al] 
l2004h is a robotic telescope located at the ESO Cerro La Silla ob- 
servatory (Chile). The REM telescope has a Ritchey-Chretien con- 
figuration with a 60-cm f/2.2 primary and an overall f/8 focal ratio 
in a fast moving alt-azimuth mount providing two stable Nasmyth 
focal stations. At one of the two foci, the telescope simultaneously 
feeds, b y means of a dichroic, two cameras: REMIR for the NIR 
(Conconi et al. 2004) and ROSS (Tosti et al. 2004) for the optical. 
The cameras both have a field of view of 10 x 10 arcmin and imag- 
ing capabilities with the usual NIR (z, J, H, and K) and Johnson- 
Cous ins VRI filters. The REM software system dCovino et al.l 
l2004h is able to manage complex observational strategies in a fully 
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autonomous way. All raw optical/NIR frames obtained with REM 
telescopes were reduced following standard procedures. Instrumen- 
tal magnitudes were obtained via aperture photometry and absolute 
calibration has been perfo rmed by means of seco ndary standard 
stars in the field reported in lHamuv fc Mazal jl989l) or by 2MAS^| 
objects in the field. The data presented here were obtained during 
2008 August-2010 April by a Guest Observer programme for an- 
nouncement of observing time AOT18 (PI: D. Impiombato) and 
AOT19 (PI: F. D'Ammando) and a long term project for AOT20 
and AOT21 (PI: F. D'Ammando). 



3.3 ATOM 

Optical observations in the Johnson R filter for this campaign 
were obtained between 2008 August and 20 10 April with th e 
0.8-m optical telescope ATOM in Namibia ( Hauser et al. I l200i) . 
ATOM is operated robotically by the High Energy Stereoscopic 
System (H.E.S.S.) collaboration and obtains automatic observa- 
tions of confirmed or potentially 7-bright blazars. Data analysis 
(debiassing, flat fieldin g, and photometry with Source-Extractor, 
iBertin & Amoutslll996l) is conducted automatically using our own 
pipel ine. For different i al pho tometry, the reference stars 2, 3, and 6 
from lHamuv fc Mazj ( Il989l) were used. 



3.4 SMA 

The 230 GHz (1.3 mm) light curve was obtained at the Submil- 
limeter Array (SMA) on Mauna Kea (Hawaii) from 2008 August 1 
to 2010 March 9. PKS 0537-441 is included in an ongoing moni- 
toring programme at the SMA to determine the fluxes of compact 
extragalactic radio sources that can be used as calibrators at mm 
wavelengths . Details of the obse rvations and data reduction can 
be found in ^ GurweU et al.l ( l2007h . Data from this programme are 
updated regularly and are available at SMA websitqj. Additional 
SMA data are from a programme led by A. Wehrlqj to monitor 
fluxes of blazars on the Fermi LAT Monitored Source LisQ. 



4 GAMMA-RAY SPECTRAL AND TEMPORAL 
VARIABILITY 

4.1 Spectral behaviour 

The 7-ray spectral analysis of PKS 0537-441 was performed both 
for the first 20-month (2008 August 4-2010 April 4; MJD 54682- 
55290) of Fermi-LAT observations, and for 5 sub-periods: period 
I (2009 February 28-May 11; MJD 54890-54962), period 2 (2009 
June 6-July 30; MJD 54988-55042), period 3 (2009 August 10- 
October 27; MJD 55053-55 13 1), period 4 (2009 November 6- 
December 9; MJD 55 141-55 174), and period 5 (2010 February 
28-April 4; MJD 55255-55290). The sub-periods, reported also 
in Fig. [T] have been chosen for selecting different activity states 
in 7 rays. The results are reported in Table [T] A hardening of the 
7-ray spectrum has been observed during period 5, which corre- 
sponds to the highest activity state. This "harder when brighter" 



* http://www.ipac.caltech.edu/2mass/ 

^ http://smal.sma.hawaii.edu/callist/callist.html. Use in publication re- 
quires obtaining permission in advance. 
^ awehrle@spacescience.org 

^ http://fermi.gsfc.nasa.gov/ssc/data/policy/LAT_Monitored_Sources.html 
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Figure 2. 20-month 7-ray spectrum of PKS 0537^41 measured by Fermi- 
LAT between 2008 August 4 and 2010 April 4. Solid brown line: power-law 
spectrum fit. Dashed orange line: broken power-law spectrum fit. Dotted 
magenta line: log-parabola spectrum fit. 



behaviour was already reported in other bright flat spectrum ra- 
dio quasars (FSRQ s) and low-synchrotron-peaked (LSP) BL Lacs 
jAbdo et I l2010dl) , even if only moderate variation (AF < 0.3) 
has been observed. By contrast we noted that during period 2 no 
hardening of the spectrum of PKS 0537—441 has been observed 
despite the high flux, indicating that the "harder when brighter" 
effect is not detected in all high activity periods of the source. The 
difference could be due to different causes or locations of the differ- 
ent flares. We noted that a "harder when brighter" trend is not a uni- 
versal behaviour in blazar 7-ray flares. For example, in 3C 454.3, 
there is only very weak evidence of a "har der when brighter" trend 
jAckermann et al .l2OI0l ;l Abdo et al I20T ij) and there is no evidence 
for signi ficant correlation betw een 7-ray flux and photon index in 
3C 279 ( iHavashida et"aLll 20121) . The 20-month spectrum is shown 
in Fig.|2] The energy spectrum was built by dividing the whole en- 
ergy range (0.1-100 GeV) into bands, requiring TS > 50 and/or 
more than 8 photons predicted by the maximum likelihood source 
for each bin except for the last one. This results in 14 energy bins 
and one upper limit for the 20-month spectrum, and 3 bins and one 
upper limit for the 5 spectra in the sub-periods. For each energy 
bin a maximum likelihood analysis, fixing the spectral index at the 
respective global likelihood values computed in the entire 0.1-100 
GeV energy range, was performed. 

By visual inspection the 20-month spectrum seems to be 
curved (Fig. |2}. To explore this curvature the spectrum was fit- 
ted in the 0.1-100 GeV energy band with two alternative spec- 
tral models with respect to the simple power-law: a broken power- 
law (BPL), dN/dE oc (B/E'o)"^! for E < Eb, dN/dE oc 
(E/Eo)^^^ for E > Et, where Fi and r2 are the photon in- 
dices below and above the energy break Et; a LogParabola (LP), 
dN/dE cx (_B/£o)"""''''°®<^/^''\ where the parameter a is the 
spectral slope at a refere nce energy Ep and f3 measures the curva - 
ture around the peak (see lLandau et aI.lll986l ; lMassaro et al.ll2004h . 
In order to estimate as accurately as possible the energy break pa- 
rameter of the BPL function, we studied the profile of the likeli- 
hood function, fixing Et at different values between 100 MeV and 
5 GeV with a step of 50 MeV. This likelihood profile was fitted with 
a parabolic function; the minimum for the parabola corresponds to 
Et = (2290 ± 390) MeV. Fixing Et at this value during the BPL 
likelihood fit, the resuks are: prefactor = (3.31 ± O.I3)x 10~^^, Fi 
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Table 1. Results of the spectral fits to the Ferm(-LAT data in the 0.1-100 
GeV energy range for the 20-month and the 5 sub-periods with a power-law 
model. 



Period 


Flux (E > 100 MeV) 


r 


TS 




[xlO-** ph s-i] 






20-month 


39.9 ± 0.8 


2.20 ± 0.01 


15661 


Period 1 


29.5 ± 1.7 


2.31 ±0.04 


899 


Period 2 


60.9 ± 0.2 


2.27 ± 0.02 


2247 


Period 3 


33.0 ± 0.4 


2.36 ± 0.01 


1195 


Period 4 


39.7 ± 3.1 


2.55 ± 0.08 


2514 


Period 5 


73.1 ± 2.6 


2.09 ± 0.03 


4500 




1.0 



1.5 2.0 

log,.(Al) [days] 



2.5 



Figure 3. First order structure function of the 3-day bin 7-ray light curve 
for the period 2008 August 4-2010 April 4. 




Figure 4. Periodogram P(f) of the 3-day bin 7-ray hght curve. No signifi- 
cant peak is present in P(f); therefore no significant power spectrum com- 
ponent and no periodicity is present. 
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= 2.10 ± 0.02, Ta = 2.69 ± 0.07, with TS = 15733 and an integral 
flux in the 0.1-100 GeV energy range of (37.6 ± 0.8) x 10"** ph 
cm~^ s~^. The LogParabola fit was performed fixing the reference 
energy Eq at 300 MeV; the normalization was found to be (4.28 
± 0.08) X 10"^°, a = 2.01 ± 0.03 and /3 = 0.09 ± 0.01 with TS 
= 15732 and an integral flux in the 0.1-100 GeV energy range of 
(38.7 ± 0.7)x 10"* ph cm~^ s"\ We used a likehhood ratio test 
to check the PL model (null hypothesis) against the BPL model 
(alternative hypoth esis). These values may be compared, following 
iNolan et^ l2012i) . by defining the significance of the curvature 
TScurve=TSBPL-TSpL=72 Corresponding to ~8.5 a. In the same 
way for the LP model we obtain TScurvc=71 corresponding to ~8.4 
a. This shows that the 20-month LAT spectrum of PKS 0537-441 
does not follow a power-law function, and the curvature is signifi- 
cant. 



Figure 5. The Discrete Autocorrelation Function applied to the 3-day bin 
7-ray light curve. DACF points out a zero-value crossing time lag between 
69 and 72 days. 



days, probably due to both the contamination of finite-series edge 
effects and to the main outburst and active phase of PKS 0537- 
441, which can be seen in the central part of the light curve from 
about MJD 54885 to about 55087. A hint of a possible moderate 
break is found at timescales between 48 and 51 days, but this does 
not necessarily imply a local characteristic time scale, considering 
that the periodogram of the light curve does not show significant 
evidence for any characteristic and recurrent time-scale (Fig. |4ll. 
The discrete autocorrelation function (DACF) does not show peaks 
at time lags below the 200-day limit, and the zero-crossing time lag 
is placed between 69 and 72 days (Fig.|5). 



4.2 Temporal behaviour 

To investigate the temporal activity of PKS 0537-441, we used the 
light curve of the integrated 7-ray flux over the entire period 2008 
August 4-2010 April 4 with 3-day time bins in the band 0.1-100 
GeV already presented in Sect 2. 

In Fig. [3] the first order structure function (SF) of the 7-ray 
light curve is reported. The SF analysis shows a power-law increas- 
ing trend from zero-lag to (147 ± 3) days lag, resulting in a 1//" 
power density spectrum (PDS) with a = 1.5 ± 0.2. This trend 
implies universality from 3-day to about 5-month time-scales, and 
the index value represents a temporal variability placed halfway 
between the flickering (red noise) and the shot noise (Brownian- 
driven process) behaviour. This slope value is in agreement with 
the average value of the PDS directly evaluated in the frequency 
domain using the 3-day bin light curves extract ed for the 9 brighte st 
FSRQs of the first year of Fermi observations ( lAbdo et al.l2010cr) . 

A consistent drop is seen in Fig.[3]at lags longer than about 200 



5 DISCRETE CROSS CORRELATION ANALYSIS 

Correlations between the best sampled multiwavelength light 
curves of PKS 0537—441 were investigated by comput- 
ing the discrete cross correlation function (DCCF), following 
lEdelson fc Kroliki fl988i) . Lag and correlation strength were com- 
puted by fitting a gaussian profile to the correlation peak in the 
DCCF and uncertainties in these parameters were estimated by a 
Monte Carlo method (see Peterson et al. 1998). DCCFs were cal- 
cul ated both wit h detrended and non-detrended light curves (see 
e.g. lWelshlll999l for details about the detrending procedure) in or- 
der to look for differences in correlation of the rapid variations and 
the slower long term variability. 

The combined ATOM and REM i?-band data provide the best 
overlap with the LAT 7-ray light curve. The DCCF between these 
two data sets, without detrending, is shown in Fig. |6] left panel 
(here and in the following positive lag in 7 rays - R means that 7- 
ray flux variations lead those in the _R-band). The peak in Fig.|6](left 
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Figure 6. Discrete Cross Con'elation Function between tlie 7-ray and i?-band (left panel), 7-ray and _ff -band (center panel), i?-band and _ff -band (right panel) 
light curves. In each case the DCCF shows a correlation peak consistent with lag = 0. 



Table 2. Time lag and strength of the cross coiTelation peak estimated by 
a gaussian fit. All data were used in the case with no detrending, while the 
detrended light curves used data from MJD 54890 to 55290. The correlation 
values have been corrected for the effect of white noise due to measurement 



2009 



2009.5 



2010 







Not detrended 


Detrended 


DCCF 




Lag (days) 


Max 


Lag (days) 


Max 


7 rays - 


R 


3.1 ± 1.9 


0.92 ± 0.06 


0.5 ± 2.4 


0.57 ± 0.07 


7 rays - 


K 


1.9 ±2.5 


0.54 ± 0.07 


0.1 ±2.6 


0.43 ± 0.08 


R-K 




0.0 ± 1.1 


0.71 ± 0.06 


-0.5 ± 0.9 


0.64 ± 0.05 



panel) corresponds to a correlation of 0.92 ± 0.06 when corrected 
for the effect of measurement errors (which reduces the correla- 
tion). A Spearman rank coiTelation test between the 7-ray and R- 
band flux gives coefficient rs of 0.6609 with a probability of chance 
occurence < 10~^, confirming the positive correlation. The corre- 
lation between the 7-ray and the K-hwA is weaker, 0.54 ± 0.07 
(Fig-ID center panel). This is not an effect of the sparser sampling 
of the K-bwiA light curve, since the difference is present also when 
the analysis is limited to the time range where both bands are well 
sampled (MJD 54890-55290). Furthermore the DCCFs show no 
significant time lag, with an uncertainty of a few days, between the 
R and K bands (Fig. [6] right panel). 

It is clear from the light curves that there is a correlation on 
long time scales and this is also shown by the results of the DC- 
CFs for the detrended light curves (see Table. |2). The correlations 
are weaker compared to the non-detrended ones and also here no 
significant time lag was detected. The position of the peak is con- 
sistent with zero time lag (1.3 ± 1.5 days). When both light curves 
were detrended with a second order polynomial the correlation was 
reduced to 0.51 ± 0.07, but still consistent with zero time lag (0.3 
± 1.8 days). However different behaviours are quite evident going 
into details of the NIR, optical, and 7-ray light curves (see Figs. |7] 
and [8j. A similar situation has been already observed for other 
brigh t 7-ray blazars (e.g. 3C 454.3, 4C +38.41: lRarteri et al.l201ll 
I2OI2I) . This complex behaviour could be in agre ement with the tur - 
bulent extreme multi-cells scenario proposed bv lMarscher I ( |201^ . 
In this context the slope of the PDS obtained for the 7-ray light 
curve indicates a significant contribution of red noise and thus of a 
random process such as the turbulence that is modulating the emis- 



6 X-RAY TO MM BEHAVIOUR 

The BAT spectrum of the data collected during 2004 November- 
2009 August was fit with a power law (dN/dE oc ^-rBAT) ^ 
photon index Tbat = 1.5 ± 0.5 and a 15-150 keV flux of (1.5lo;5) 
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Figure 7. Comparison between 7-ray and ii'-band light curves. We super- 
impose 7-ray (black triangles) and iC-band (red empty circles) light curves 
normalizing 7-ray and K flux values over the whole observing period to the 
respective peak flux values. The time binning is 3 days for 7-ray data and 
1-day for iC-band data set. 



X 10"^^ erg cm"^ s"^ (x?od = 0.968, with 6 d.o.f.). As a compar- 
ison the same photon index with a flux a factor of two higher was 
obtained for the BAT spectrum extracted from 2004 December to 
2005 November, dominated by episodes of more intense activity 
jPian et alJboOTi) . indicating a quite stable spectral shape in hard 
X-rays. 

During 2008 October - 2010 March SwiftlXKl observed the 
source with a 0.3-10 keV flux in the range (3.2-7.7) x 10^^^ erg 
cm~^ s~^, with a photon index varying in the range 1.6-1.9. The 
interval of photon index values is very similar to that previously 
observed by SwiftlXKI spanning a larger flux range (Fx = 1.6- 
1.8, jPian^t alj|2007h . A Spearman rank correlation test between 
the 7-ray and the X-ray flux gives a coefficient r^ of 0.6176 with a 
probability of chance occurrence of 8 x 10""', indicating a positive 
correlation. The small variability amplitude (calculated as the ratio 
of maximum to minimum flux) observed in X-rays (~ 2.5) with 
respect to the 7 rays (~ 11) could be an indication that the X-ray 
emission is produced by the low-energy tail of the same electron 
distribution. The peak flux observed in 2008 October is a factor of 
two lower with respect to the 2004-2005 SwiftKKI observations 
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Figure 8. Compaiison between 7-ray and ij-band light curves. We super- 
impose 7-ray (black triangles) and ij-band (red empty circles) REM light 
curves normalizing 7-ray and R flux values over the whole observing pe- 
riod to the respective peak flux values. The time binning is 3-day for 7-ray 
data and 1-day for ij-band data set. 
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Figure 9. Swift/XRT photon index of PKS 0537-441 as a function of the 
0.3-10 keV flux. The observations performed on 2008 October are shown 
with red circles. 



jPian et alJboOTi) . Interestingly the high X-ray flux in 2008 Octo- 
ber coincides with a low activity state both in NIR-to-UV and 7-ray 
bands, suggesting that a second emission component can give a sig- 
nificant contribution in soft X-rays in that period (e.g. an increase 
of the synchrotron self-Compton component). We noted that X-ray 
flux variations had n o optical/UV counterpart also in 2005 Novem- 
ber jPian et al.|[2007D . Similarly an X-ray flare w ith no counterpart 
at other frequencies ha s been observed in 3C 279 dAbdo et al.l2"oTol : 
iHavashida et alj2012l) . 

Fig. |9] shows the XRT photon indices as a function of the 
fluxes in the 0.3-10 keV band; no significant change in the X-ray 
spectrum with the increase in the flux has been observed, and thus 
there is no hint of an hard-spectrum additional component becom- 
ing persistent and significant at higher fluxes. However we note 
a relatively hard spectrum (F ~1.6-1.7) during the 2008 October 
observations (red circles in Fig. |9). The high correlatio n between 
the X- ray and optical (V^-band) hght curves reported in lPian et al.l 
is not completely confirmed by our long-term monitoring of 
the source. In fact the peak of the activity in l/-band observed by 
SwiftfUYCfT in 2010 March has no counterpart in the X-ray light 
curve. This, together with the increase of the flux observed in X- 
rays in October 2008 with no counterpart at other wavelengths, 
suggests a more complex connection between optical. X-ray (and 7 
ray) emission over a long period. However the sparse sampling of 
the X-ray light curve does not allow us to make a precise compari- 
son with the NIR, optical as well as 7-ray light curves. By contrast, 
the ii'-band and i?-band are well sampled and a comparison with 
the 7-ray light curve showed an overall correlation (see Sect. 5). 
However on inspecting the light curves on short time scales we 
note a significant increase of the NIR and optical fluxes in 2009 
September and November with a smaller variability amplitude in 
7 rays (see Figs. |7] and [§}, peaking on 2009 September 21 (MJD 
55095) and November 12 (MJD 55147) (see Fig.[D. The largest 



amplit ude variations are usually detected in 7 rays dAbdo et al.l 
l2010d) , and thus this behaviour is quite peculiar. Recently, a dif- 
ferent behaviour between 7-ray and NIR/optical light curves has 
been observed in the intermediate-frequency-peaked BL Lac 3C 
66 A in 2008 and 2009-2010, with a good correlation in 2008 and 
an increasin g NIR/optical flu x with no counterpart in 7 rays during 
2009-2010 dltoh et aljlloi^ . In the same way, notwithstanding a 
significant correlation estimated between optical and 7-ray bands 
during 2008-2011 a strong optical flare without a 7-ray counter- 
part has been dete cted from the FSRQ 4C +38.41 in 2011 July 
dRaiteri et all2012[) . 

During 2008-2010 an increase with an amplitude of a factor 
of ~7 was observed in IZ-band, significantly lower with respect 
to 2004-2005 observations (a factor of ~60), confir ming the huge 
flaring activity observ ed by REM and Swift in 2005 dDolcini et alj 
l2005l : 1pran et al.l2007h . A decreasing variability amplitude was ob- 
served from NIR (a factor of ~12-13) to UV (a factor of ~5). In 
particular a clear change of variability amplitude was observed be- 
tween 7?-band (a factor of ^--^14) and l/-band (a factor of --^7). 

There appears to be a good correlation between the 230 GHz 
light curve collected by SMA and the 7-ray light curve, in particu- 
lar a similar increase was observed in 2010 reaching a peak value of 
8.05 Jy on 2010 March 9 for SMA during the brightest 7-ray flaring 
period. The lack of SMA observations after this date does not allow 
us to determine if the 230 GHz flux density is still increasing and 
whether the mm-peak is strictly simultaneous with the 7-ray peak. 
However the contemporaneous increase in the two energy bands in 
2010 seems to indicate that the emission region at mm and 7 rays is 
co-spatial and thus the 7-ray flaring activity probably originates at 
large distance from the central engine, beyond the broad line region 
(BLR). 
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Figure 10. Modeling of the SEDs of PKS 0537-441 in a low, average, and 
high state including Fermi, Swift, ATOM, REM, and SMA data. See text for 
details. 



7 MODELING THE SEDS 

We build three SEDs for PKS 0537-441 in three different activity 
states: low, average, and high. We used REM, ATOM, and Swift 
data collected on 2009 March 4 (MJD 54894), 2009 June 4 (MJD 
54986), and 2010 March (MJD 55258), and LAT spectra calculated 
over the period 2009 February 28 - May 1 1 (MJD 54890-54962), 

2008 August 4 - 2010 February 4 (MJD 54682-55231), and 2010 
February 28 - April 4 (MJD 55255-55290) for the low, average, 
and high activity states, respectively. For each SED we used the 
SMA observation available nearest to the Swift observation, i.e. 

2009 February 27, 2009 July 15, and 2010 March 7. For the high ac- 
tivity state also the Planck and Wide-Field Infrared S urvey Explorer 
(WISE ) data collected on 2010 March 7-22 (from iGiommi et alj 
l2012h and March 4-7, respectively, has been reported. Finally we 
report in the SED the average BAT spectrum. 

All three SEDs indicate a lower-frequency component peaking 
at ~ 10^^ — 10^* Hz and a higher-frequency component peaking 
around 10^^ - 10^^ Hz. We attempted to reproduce the SEDs 
with leptonic models which in clude synchrotron and synchrotron 
self-Compton (SSC) scattering jFinke et al.ll2008h . These attempts 
failed, as expected for FSRQs and some LSP BL Lacs. We then 
added an external Compton (EC) component of seed photons from 
a dust torus. We were able to obtain reasonable fits to all three 
states, with the only difference between states being the electron 
distribution. These model fits are shown as curves in Figure [TO] and 
the m odel parameters can be found in Table |3] See .Dermer et al.i 
( l2009h for a description of the model parameters. The dust torus 
was modelled as a one-dimensional ring around the black hole, 
aligned orthogonal to the jet, and its parameters are approximately 
in agreement wit h dust emitting at the sublimation radius with the 
formula given bv lNenkova et al.l ( 1200 8h . The dust luminosity in the 
model is rather low, but note that there is no hint of dust or disc 
emission in the IR-optical por tion of the SED, in agreement with 
the results reported in Impiomb ato et al.l (201 1 ). For all states, we 
found that an electron distribution with two power laws (a single 
break) was not enough to explain the observed SED; an additional 
power law and break were needed, in particular for reproducing the 
NIR-optical part of the spectrum, so that the electron distribution is 



given by 



Tmin ^ T ^ 'ybrk,l 
; I'brk,! < 7' < l'brk,2 



« 7 , 7brfc,2 < 7 < 7rT 



For the modeling, the variability time was chosen to be a bit more 
than two days, consistent with the observed light curve. This con- 
strains the size of the emitting region, given a Doppler factor. In 
our model fits, the primary emitting region is optically thin to syn- 
chrotron self-absorption down to ~ 200 GHz, so that it can re- 
produce the data at higher frequencies. Below this frequency, pre- 
sumably th e radio emissi on comes from other, larger emitting re- 
gions (e.g., lKonigllll98ll) . In our model fits, this blob is quite far 
from the black hole, J5 1 pc. If the emitting blob takes up the en- 
tire cross section of the jet, then it would have an opening angle 
of Sopen ~ 2° . This is approximately consistent with the opening 
angles f ound for other blaz ars based on multi-epoch VLBI obser- 
vations jjorstad et al]|2005l) . 

In the "average state", the LAT spectrum shows a clear devia- 
tion from a single power-law (see Section 3). This has been foun d 
for numerous other FSRQs and LSP BL Lacs jAbdo et al.ll2010dh . 
most notably for the extremely 7-ray bright FSRQ 3C 454.3 
jAbdo et al.l(2009l : lAckermann et al.ll201ol; lAbdo et al.ll201 ih . The 
cause of this curvature, often characterized as a spectral break 
(when fit with a double power-law) is a bit of a mystery. It has 
been suggested that it is due to a feature in the electron distribu- 
tion jAbdo et alj|200^ : due to 77 abs orption of 7-ra ys with He II 
Lyg b road line photons jPoutanen & S tem 2010; Stem & PoutanenI 
l201lh ; due to a combination of Compton scattering of two seed 
photon sources, for example, d irectly from the accret ion disc, 
and from broad line emission jFinke & Denned I2OIOI) ; or due 
to Klein-Nishin a effects from C ompton-scattering of H I Lya 
photons ( Acken nann et al.l I2OIO ). In the "high state" SED for 
PKS 0537—441, there is clearly a spectral break in the IR-optical 
spectmm, around 3 x 10^* Hz (~ 1 /im). The combination of a 
break found in both the synchrotron and Compton-scattered spec- 
tmm seems to indicate a common origin for the break, most likely 
due to an intrinsic feature in the underlying electron distribution. 
Indeed, that is how we model it for this source. Although the broken 
IR/optical spectmm in the "high state", and the curved 7-ray spec- 
tmm in the "average state" are not contemporaneous, there does 
seem to be a hint of curvature in the "high state" LAT spectmm, 
especially when the upper limit in the highest energy bin is con- 
sidered. It is certainly possible, maybe even likely, that the cause 
of the 7-ray spectral breaks in all blazars which show them is re- 
lated to features in the electron distribution, although one should 
be careful not to over-generalise. Simultaneous IR, optical and 7- 
ray observations of other blazars could be very useful for mak- 
ing this det ermination, although c ontamination in the IR by a dust 
toms (e.g., 'Malmrose etal."201 F) or another synchrotron compo- 
nent (e.g.,fOgle et al. 201 1) could make this difficult. 

The three states differ very little in the radio and X-rays. This 
emission is caused by synchrotron and EC emission, respectively, 
of the lower portion of the electron distribution, although there is 
some additional SSC emission for the X-rays. Consequently, the jet 
power in electrons (see Table |3) varies little between states, since 
the pi = 2.0 portion of the electron distribution contains most 
of the electrons' energy. The model implies emission very close to 
equipartition between electrons and Poynting flux for all model fits, 
with a slight dominance of Poynting flux. 
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Table 3. Model parameters for the SEDs of PKS 0537-441. 



Parameter Symbol Low State Average State High State 



Redshift 


z 






0.896 


0.896 


0.896 


Bulk Lorentz Factor 


r 






50 


50 


50 


Doppler Factor 








50 


50 


50 


Magnetic Field 


B 






0.2 


0.2 


0.2 G 


Variability Time-Scale 






2 X 10^ 


2 X 10^5 


2 X 10^ 


Comoving radius of blob 




1, 


.6 


X 10^'' cm 


1.6 X IQi'^ cm 


1.6 X lO^'^ cm 


Jet Height 


r 


3, 


.7 


X lO'^^ cm 


3.7 X IQi* cm 


3.7 X 10^8 cm 


First Electron Spectral Index 


Pi 






2.0 


2.0 


2.0 


Second Electron Spectral Index 


P2 






3.1 


3.1 


3.1 


Third Electron Spectral Index 


P3 






3.9 


4.6 


4.6 


Minimum Electron Lorentz Factor 


7' ■ 

' jmn 






4.0 


1.0 


3.0 


Break Electron Lorentz Factor 


'ybrk,! 




3 


.0 X 10^ 


4.7 X 10^ 


7.0 X 10^ 


Break Electron Lorentz Factor 






9 


.0 X 10^ 


2.6 X 10=* 


4.2 X 10^ 


Maximum Electron Lorentz Factor 


'^max 




1 


.3 X 10"* 


2.3 X IQ-* 


3.0 X 10"' 


Dust Torus luminosity 


^dust 


3.3 


X 


lO'*^ erg s-i 


3.3 X 10*2 ergs-i 


3.3 X 10''2 ergs-l 


Dust Torus temperature 


'-^dust 




3 


X 10^ K 


3 X 10^ K 


3 X 10^ K 


Dust Torus radius 


^dust 


9 


.4 


X 10^* cm 


9.4 X IQi* cm 


9.4 X IQi* cm 


Jet Power in Magnetic Field 




1.9 


X 


lO-'S erg s-l 


1.9 X 10^"^ ergs-l 


1.9 X 10''^ ergs-l 


Jet Power in Electrons 


p. 


1.2 


X 


lO-'S erg s-l 


1.2 X 10-*^ ergs-l 


1.3 X 10''^ ergs-i 



8 CONCLUSIONS 

We presented multiwavelength observations of PKS 0537—441 
during a period of 20 months (2008 August~2010 April) including 
Fermi, Swift, ATOM, REM, and SMA data. Strong variability has 
been observed in 7 rays, with two major flaring episodes (2009 July 
and 2010 March) and a harder-when-brighter spectral behaviour, 
quite common for FSRQs and LSP BL Lacs, in 2010 March. The 
average LAT spectrum accumulated over 20 months showed a sig- 
nificant curvature, well described both by a BPL model with Et = 
(2290 ± 390) MeV and a log-parabola model. A quite steep PDS 
slope (a = 1.5 ± 0.2) has been es timated, in agreement with those 
observed for the brightest FSRQs dAbdo et al. 2010c), suggesting a 
random-walk underlying mechanism. 

Only moderate variability has been observed in X-rays with 
no correlation between flux and photon index. An increase of the 
flux in the X-ray band with no counterpart at the other wavelengths 
has been observed in 2008 October, suggesting a significant con- 
tribution of a second component (e.g. SSC emission) in some high 
activity states. No clear correlation of the X-ray emission with the 
NIR or optical bands has been found. Only a correlation with the 
gamma-ray band has been observed. In addition we observed cor- 
relation between the 7-ray band with the ii-band and more weakly 
with the ii'-band, with no significant time lag. However, together 
with this "overall" correlation we note that on finer timescales the 
light curves showed differences. In particular, two flaring episodes 
in NIR and optical bands have been observed in 2009 Septem- 
ber and even more strongly in 2009 November with no significant 
counterparts in 7 rays. 

As seen for FSRQs and LSP BL Lacs, the SED of the source 
cannot be modelled by a simple SSC model. We included an EC 
component of seed photons from a dust torus. The 230 GHz light 
curve showed an increase simultaneously with the 7 rays, suggest- 
ing co-spatiality of the mm and 7-ray emission regions at large dis- 
tance from the central engine and thus the dust torus as the possible 
main source of seed photons. The low, average, and high activity 
SEDs of the source could be fit by changing only the electron dis- 
tribution parameters. This has been observed previously for Mrk 



501 jPetrv et al.l l2000). although this is a high-synchrotron-peaked 
BL Lac, while PKS 0537-441 is an LSP BL Lac/FSRQ. In our mod- 
eling, we have found that two breaks in the electron distribution are 
necessary. A spectral break in the NIR-optical spectrum has been 
found, in agreement also with a significant change of the flux vari- 
ability amplitude below and above the _R-band. This break together 
with the curvature observed in 7 rays (possibly characterized as a 
spectral break) seems to indicate a common origin for them, most 
likely due to an intrinsic feature in the underlying electron distribu- 
tion. 

Broad emission lines hav e been observed during a low ac- 
tivity state of PKS 0537-441 IPian et al.ll2005l) . in contrast with 
the initial classification as a BL Lac object. However, recent stud- 
ies of a large sample of blazars showed that a classification de- 
pending only on their optical/UV spectral pro perties is not effi- 
cient to distinguish FS RQs and BL Lac objects dAbdo et aPbOlOal: 
lGhiselliniet"al]l201lh . Variability and spectral properties in 7 rays 
indicate a FSRQ-like behaviour, in agreement also with the SED 
properties. The observed isotropic 7-ray luminosity in the 0.1-100 
GeV energy range is 8.2 x 10*^ erg s^^ for the 20-month interval 
considered, reaching a peak value of 2.6x10''* erg s~' on 3-day 
time scale at the end of 2010 March. These values are comparable 
to those of the bright FSRQs jGhisellini et alj2O10l) . 

The multifrequency observations presented here give new 
clues, but also offer new questions on the astrophysical mecha- 
nisms at work in PKS 0537—441. Further mm to 7-ray observa- 
tions will be fundamental to investigate in even more detail the 
correlations by using the long-term variability in different energy 
bands to achieve a complete understanding of the fundamental pro- 
cesses that underlie the behaviour of this source. 
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